Genome-wide bioinformatics analysis of human protease capacity for proteolytic cleavage of the SARS-CoV-2 spike glycoprotein.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) primarily enters the cell by binding the virus's spike (S) glycoprotein to the angiotensin-converting enzyme 2 receptor on the cell surface, followed by proteolytic cleavage by host proteases. Studies have identified furin and transmembrane protease serine 2 proteases in priming and triggering cleavages of the S glycoprotein, converting it into a fusion-competent form and initiating membrane fusion, respectively. Alternatively, SARS-CoV-2 can enter the cell through the endocytic pathway, where activation is triggered by lysosomal cathepsin L. However, other proteases are also suspected to be involved in both entry routes. In this study, we conducted a genome-wide bioinformatics analysis to explore the capacity of human proteases in hydrolyzing peptide bonds of the S glycoprotein. Predictive models of sequence specificity for 169 human proteases were constructed and applied to the S glycoprotein together with the method for predicting structural susceptibility to proteolysis of protein regions. After validating our approach on extensively studied S2' and S1/S2 cleavage sites, we applied our method to each peptide bond of the S glycoprotein across all 169 proteases. Our results indicate that various members of the proprotein convertase subtilisin/kexin type, type II transmembrane family serine protease, and kallikrein families, as well as specific coagulation factors, are capable of cleaving S2' or S1/S2 sites. We have also identified a potential cleavage site of cathepsin L at the K790 position within the S2' loop. Structural analysis suggests that cleavage of this site induces conformational changes similar to the cleavage at the R815 (S2') position, leading to the exposure of the fusion peptide and subsequent fusion with the membrane. Other potential cleavage sites and the influence of mutations in common SARS-CoV-2 variants on proteolytic efficiency are discussed.IMPORTANCEThe entry of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) into the cell, activated by host proteases, is considerably more complex in coronaviruses than in most other viruses and is not fully understood. There is evidence that other proteases beyond the known furin and transmembrane protease serine 2 can activate the spike protein. Another example of uncertainty is the cleavage site for the alternative endocytic route of SARS-CoV-2 entrance, which is still unknown. Bioinformatics methods, modeling protease specificity and estimating the structural susceptibility of protein regions to proteolysis, can aid in studying this topic by predicting the involved proteases and their cleavage sites, thereby substantially reducing the amount of experimental work. Elucidating the mechanisms of spike protein activation is crucial for preventing possible future coronavirus pandemics and developing antiviral drugs.

Bioactive glycans in a microbiome-directed food for children with malnutrition.

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1-4. Here we analyse biospecimens from a randomized, controlled trial of a microbiome-directed complementary food (MDCF-2) that produced superior rates of weight gain compared with a calorically more dense conventional ready-to-use supplementary food in 12-18-month-old Bangladeshi children with moderate acute malnutrition4. We reconstructed 1,000 bacterial genomes (metagenome-assembled genomes (MAGs)) from the faecal microbiomes of trial participants, identified 75 MAGs of which the abundances were positively associated with ponderal growth (change in weight-for-length Z score (WLZ)), characterized changes in MAG gene expression as a function of treatment type and WLZ response, and quantified carbohydrate structures in MDCF-2 and faeces. The results reveal that two Prevotella copri MAGs that are positively associated with WLZ are the principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilizing the component glycans of MDCF-2. The predicted specificities of carbohydrate-active enzymes expressed by their polysaccharide-utilization loci are correlated with (1) the in vitro growth of Bangladeshi P. copri strains, possessing varying degrees of polysaccharide-utilization loci and genomic conservation with these MAGs, in defined medium containing different purified glycans representative of those in MDCF-2, and (2) the levels of faecal carbohydrate structures in the trial participants. These associations suggest that identifying bioactive glycan structures in MDCFs metabolized by growth-associated bacterial taxa will help to guide recommendations about their use in children with acute malnutrition and enable the development of additional formulations.

Bioactive glycans in a microbiome-directed food for malnourished children.

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1-4. Designing effective microbiome-directed therapeutic foods to repair these perturbations requires knowledge about how food components interact with the microbiome to alter its expressed functions. Here we use biospecimens from a randomized, controlled trial of a microbiome-directed complementary food prototype (MDCF-2) that produced superior rates of weight gain compared to a conventional ready-to-use supplementary food (RUSF) in 12-18-month-old Bangladeshi children with moderate acute malnutrition (MAM)4. We reconstructed 1000 bacterial genomes (metagenome-assembled genomes, MAGs) present in their fecal microbiomes, identified 75 whose abundances were positively associated with weight gain (change in weight-for-length Z score, WLZ), characterized gene expression changes in these MAGs as a function of treatment type and WLZ response, and used mass spectrometry to quantify carbohydrate structures in MDCF-2 and feces. The results reveal treatment-induced changes in expression of carbohydrate metabolic pathways in WLZ-associated MAGs. Comparing participants consuming MDCF-2 versus RUSF, and MDCF-2-treated children in the upper versus lower quartiles of WLZ responses revealed that two Prevotella copri MAGs positively associated with WLZ were principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilization of its component glycans. Moreover, the predicted specificities of carbohydrate active enzymes expressed by polysaccharide utilization loci (PULs) in these two MAGs correlate with the (i) in vitro growth of Bangladeshi P. copri strains, possessing differing degrees of PUL and overall genomic content similarity to these MAGs, cultured in defined medium containing different purified glycans representative of those in MDCF-2, and (ii) levels of carbohydrate structures identified in feces from clinical trial participants. In the accompanying paper5, we use a gnotobiotic mouse model colonized with age- and WLZ-associated bacterial taxa cultured from this study population, and fed diets resembling those consumed by study participants, to directly test the relationship between P. copri, MDCF-2 glycan metabolism, host ponderal growth responses, and intestinal gene expression and metabolism. The ability to identify bioactive glycan structures in MDCFs that are metabolized by growth-associated bacterial taxa will help guide recommendations about use of this MDCF for children with acute malnutrition representing different geographic locales and ages, as well as enable development of bioequivalent, or more efficacious, formulations composed of culturally acceptable and affordable ingredients.

Predicting Structural Susceptibility of Proteins to Proteolytic Processing.

The importance of 3D protein structure in proteolytic processing is well known. However, despite the plethora of existing methods for predicting proteolytic sites, only a few of them utilize the structural features of potential substrates as predictors. Moreover, to our knowledge, there is currently no method available for predicting the structural susceptibility of protein regions to proteolysis. We developed such a method using data from CutDB, a database that contains experimentally verified proteolytic events. For prediction, we utilized structural features that have been shown to influence proteolysis in earlier studies, such as solvent accessibility, secondary structure, and temperature factor. Additionally, we introduced new structural features, including length of protruded loops and flexibility of protein termini. To maximize the prediction quality of the method, we carefully curated the training set, selected an appropriate machine learning method, and sampled negative examples to determine the optimal positive-to-negative class size ratio. We demonstrated that combining our method with models of protease primary specificity can outperform existing bioinformatics methods for the prediction of proteolytic sites. We also discussed the possibility of utilizing this method for bioinformatics prediction of other post-translational modifications.

Glycidamide-induced hypermutation in yeast single-stranded DNA reveals a ubiquitous clock-like mutational motif in humans.

Mutagens often prefer specific nucleotides or oligonucleotide motifs that can be revealed by studying the hypermutation spectra in single-stranded (ss) DNA. We utilized a yeast model to explore mutagenesis by glycidamide, a simple epoxide formed endogenously in humans from the environmental toxicant acrylamide. Glycidamide caused ssDNA hypermutation in yeast predominantly in cytosines and adenines. The most frequent mutations in adenines occurred in the nAt→nGt trinucleotide motif. Base substitutions A→G in this motif relied on Rev1 translesion polymerase activity. Inactivating Rev1 did not alter the nAt trinucleotide preference, suggesting it may be an intrinsic specificity of the chemical reaction between glycidamide and adenine in the ssDNA. We found this mutational motif enriched in published sequencing data from glycidamide-treated mouse cells and ubiquitous in human cancers. In cancers, this motif was positively correlated with the single base substitution (SBS) smoking-associated SBS4 signature, with the clock-like signatures SBS1, SBS5, and was strongly correlated with smoking history and with age of tumor donors. Clock-like feature of the motif was also revealed in cells of human skin and brain. Given its pervasiveness, we propose that this mutational motif reflects mutagenic lesions to adenines in ssDNA from a potentially broad range of endogenous and exogenous agents.

Molecular Heterogeneity of Pediatric AML with Atypical Promyelocytes Accumulation in Children-A Single Center Experience.

Acute promyelocytic leukemia (APL) pathogenesis is based on RARA gene translocations, which are of high importance in the diagnosis of and proper therapy selection for APL. However, in some cases acute myeloid leukemia (AML) demonstrates APL-like morphological features such as atypical promyelocytes accumulation. This type of AML is characterized by the involvement of other RAR family members or completely different genes. In the present study, we used conventional karyotyping, FISH and high-throughput sequencing in a group of 271 de novo AML with atypical promyelocytes accumulation. Of those, 255 cases were shown to carry a typical chromosomal translocation t(15;17)(q24;q21) with PML::RARA chimeric gene formation (94.1%). Other RARA-positive cases exhibited cryptic PML::RARA fusion without t(15;17)(q24;q21) (1.8%, n = 5) and variant t(5;17)(q35;q21) translocation with NPM1::RARA chimeric gene formation (1.5%, n = 4). However, 7 RARA-negative AMLs with atypical promyelocytes accumulation were also discovered. These cases exhibited TBL1XR1::RARB and KMT2A::SEPT6 fusions as well as mutations, e.g., NPM1 insertion and non-recurrent chromosomal aberrations. Our findings demonstrate the genetic diversity of AML with APL-like morphological features, which is of high importance for successful therapy implementation.

APOBEC mutagenesis is low in most types of non-B DNA structures.

While somatic mutations are known to be enriched in genome regions with non-canonical DNA secondary structure, the impact of particular mutagens still needs to be elucidated. Here, we demonstrate that in human cancers, the APOBEC mutagenesis is not enriched in direct repeats, mirror repeats, short tandem repeats, and G-quadruplexes, and even decreased below its level in B-DNA for cancer samples with very high APOBEC activity. In contrast, we observe that the APOBEC-induced mutational density is positively associated with APOBEC activity in inverted repeats (cruciform structures), where the impact of cytosine at the 3'-end of the hairpin loop is substantial. Surprisingly, the APOBEC-signature mutation density per TC motif in the single-stranded DNA of a G-quadruplex (G4) is lower than in the four-stranded part of G4 and in B-DNA. The APOBEC mutagenesis, as well as the UV-mutagenesis in melanoma samples, are absent in Z-DNA regions, owing to the depletion of their mutational signature motifs.

Lineage Conversion in Pediatric B-Cell Precursor Acute Leukemia under Blinatumomab Therapy.

We report incidence and deep molecular characteristics of lineage switch in 182 pediatric patients affected by B-cell precursor acute lymphoblastic leukemia (BCP-ALL), who were treated with blinatumomab. We documented six cases of lineage switch that occurred after or during blinatumomab exposure. Therefore, lineage conversion was found in 17.4% of all resistance cases (4/27) and 3.2% of relapses (2/63). Half of patients switched completely from BCP-ALL to CD19-negative acute myeloid leukemia, others retained CD19-positive B-blasts and acquired an additional CD19-negative blast population: myeloid or unclassifiable. Five patients had KMT2A gene rearrangements; one had TCF3::ZNF384 translocation. The presented cases showed consistency of gene rearrangements and fusion transcripts across initially diagnosed leukemia and lineage switch. In two of six patients, the clonal architecture assessed by IG/TR gene rearrangements was stable, while in others, loss of clones or gain of new clones was noted. KMT2A-r patients demonstrated very few additional mutations, while in the TCF3::ZNF384 case, lineage switch was accompanied by a large set of additional mutations. The immunophenotype of an existing leukemia sometimes changes via different mechanisms and with different additional molecular changes. Careful investigation of all BM compartments together with all molecular -minimal residual disease studies can lead to reliable identification of lineage switch.

Acute myeloid leukemia with t(X;6)9p11;q23);MYB-GATA1 and female sex: GATA1 insufficiency may be insufficient for pathogenesis.

Pediatric acute myeloid leukemia (AML) is genetically heterogenous (Olsson et al., 2016). t(X;6)(p11;q23) is a rare but recurrent chromosomal translocation in infant AML thought to be associated with male sex and basophilic differentiation (Dastugue et al., 1997). Here we report molecular characterization of AML with t(X;6)(p11;q23);MYB-GATA1 in two female infants and demonstrate preserved GATA1 expression in the sample tested. These findings further debunk a concept that this fusion was restricted to males, in whom it disrupts the only copy of the X-linked GATA1 gene, causing presumable complete loss of GATA1 function. Our data also demonstrate the power and efficiency of RNA sequencing for subclassification of leukemia on a clinically relevant timeline.

Analysis of gene expression and mutation data points on contribution of transcription to the mutagenesis by APOBEC enzymes.

Since the discovery of the role of the APOBEC enzymes in human cancers, the mechanisms of this type of mutagenesis remain little understood. Theoretically, targeting of single-stranded DNA by the APOBEC enzymes could occur during cellular processes leading to the unwinding of DNA double-stranded structure. Some evidence points to the importance of replication in the APOBEC mutagenesis, while the role of transcription is still underexplored. Here, we analyzed gene expression and whole genome sequencing data from five types of human cancers with substantial APOBEC activity to estimate the involvement of transcription in the APOBEC mutagenesis and compare its impact with that of replication. Using the TCN motif as the mutation signature of the APOBEC enzymes, we observed a correlation of active APOBEC mutagenesis with gene expression, confirmed the increase of APOBEC-induced mutations in early-replicating regions and estimated the relative impact of transcription and replication on the APOBEC mutagenesis. We also found that the known effect of higher density of APOBEC-induced mutations on the lagging strand was highest in middle-replicating regions and observed higher APOBEC mutation density on the sense strand, the latter bias positively correlated with the gene expression level.

Predictive models of protease specificity based on quantitative protease-activity profiling data.

Bioinformatics-based prediction of protease substrates can help to elucidate regulatory proteolytic pathways that control a broad range of biological processes such as apoptosis and blood coagulation. The majority of published predictive models are position weight matrices (PWM) reflecting specificity of proteases toward target sequence. These models are typically derived from experimental data on positions of hydrolyzed peptide bonds and show a reasonable predictive power. New emerging techniques that not only register the cleavage position but also measure catalytic efficiency of proteolysis are expected to improve the quality of predictions or at least substantially reduce the number of tested substrates required for confident predictions. The main goal of this study was to develop new prediction models based on such data and to estimate the performance of the constructed models. We used data on catalytic efficiency of proteolysis measured for eight major human matrix metalloproteinases to construct predictive models of protease specificity using a variety of regression analysis techniques. The obtained results suggest that efficiency-based (quantitative) models show a comparable performance with conventional PWM-based algorithms, while less training data are required. The derived list of candidate cleavage sites in human secreted proteins may serve as a starting point for experimental analysis.

Comparative immunogenicity and structural analysis of epitopes of different bacterial L-asparaginases.

BACKGROUND: E.coli type II L-asparaginase is widely used for treatment of acute lymphoblastic leukemia. However, serious side effects such as allergic or hypersensitivity reactions are common for L-asparaginase treatment. Methods for minimizing immune response on L-asparaginase treatment in human include bioengeneering of less immunogenic version of the enzyme or utilizing the homologous enzymes of different origin. To rationalize these approaches we compared immunogenicity of L-asparaginases from five bacterial organisms and performed sequence-structure analysis of the presumable epitope regions. METHODS: IgG and IgM immune response in C57B16 mice after immunization with Wollinella succinogenes type II (WsA), Yersinia pseudotuberculosis type II (YpA), Erwinia carotovora type II (EwA), and Rhodospirillum rubrum type I (RrA) and Escherichia coli type II (EcA) L-asparaginases was evaluated using standard ELISA method. The comparative bioinformatics analysis of structure and sequence of the bacterial L-asparaginases presumable epitope regions was performed. RESULTS: We showed different immunogenic properties of five studied L-asparaginases and confirmed the possibility of replacement of EcA with L-asparaginase from different origin as a second-line treatment. Studied L-asparaginases might be placed in the following order based on the immunogenicity level: YpA > RrA, WsA ≥ EwA > EcA. Most significant cross-immunogenicity was shown between EcA and YpA. We propose that a long N-terminus of YpA enzyme enriched with charged aminoacids and tryptophan could be a reason of higher immunogenicity of YpA in comparison with other considered enzymes. Although the recognized structural and sequence differences in putative epitope regions among five considered L-asparaginases does not fully explain experimental observation of the immunogenicity of the enzymes, the performed analysis set the foundation for further research in this direction. CONCLUSIONS: The performed studies showed different immunogenic properties of L-asparaginases and confirmed the possibility of replacement of EcA with L-asparaginase from different origin. The preferable enzymes for the second line treatment are WsA, RrA, or EwA.

APOBEC-Induced Cancer Mutations Are Uniquely Enriched in Early-Replicating, Gene-Dense, and Active Chromatin Regions.

An antiviral component of the human innate immune system-the APOBEC cytidine deaminases-was recently identified as a prominent source of mutations in cancers. Here, we investigated the distribution of APOBEC-induced mutations across the genomes of 119 breast and 24 lung cancer samples. While the rate of most mutations is known to be elevated in late-replicating regions that are characterized by reduced chromatin accessibility and low gene density, we observed a marked enrichment of APOBEC mutations in early-replicating regions. This unusual mutagenesis profile may be associated with a higher propensity to form single-strand DNA substrates for APOBEC enzymes in early-replicating regions and should be accounted for in statistical analyses of cancer genome mutation catalogs aimed at understanding the mechanisms of carcinogenesis as well as highlighting genes that are significantly mutated in cancer.

Horizontal gene transfer and genome evolution in Methanosarcina.

BACKGROUND: Genomes of Methanosarcina spp. are among the largest archaeal genomes. One suggested reason for that is massive horizontal gene transfer (HGT) from bacteria. Genes of bacterial origin may be involved in the central metabolism and solute transport, in particular sugar synthesis, sulfur metabolism, phosphate metabolism, DNA repair, transport of small molecules etc. Horizontally transferred (HT) genes are considered to play the key role in the ability of Methanosarcina spp. to inhabit diverse environments. At the moment, genomes of three Methanosarcina spp. have been sequenced, and while these genomes vary in length and number of protein-coding genes, they all have been shown to accumulate HT genes. However, previous estimates had been made when fewer archaeal genomes were known. Moreover, several Methanosarcinaceae genomes from other genera have been sequenced recently. Here, we revise the census of genes of bacterial origin in Methanosarcinaceae. RESULTS: About 5% of Methanosarcina genes have been shown to be horizontally transferred from various bacterial groups to the last common ancestor either of Methanosarcinaceae, or Methanosarcina, or later in the evolution. Simulation of the composition of the NCBI protein non-redundant database for different years demonstrates that the estimates of the HGT rate have decreased drastically since 2002, the year of publication of the first Methanosarcina genome. The phylogenetic distribution of HT gene donors is non-uniform. Most HT genes were transferred from Firmicutes and Proteobacteria, while no HGT events from Actinobacteria to the common ancestor of Methanosarcinaceae were found. About 50% of HT genes are involved in metabolism. Horizontal transfer of transcription factors is not common, while 46% of horizontally transferred genes have demonstrated differential expression in a variety of conditions. HGT of complete operons is relatively infrequent and half of HT genes do not belong to operons. CONCLUSIONS: While genes of bacterial origin are still more frequent in Methanosarcinaceae than in other Archaea, most HGT events described earlier as Methanosarcina-specific seem to have occurred before the divergence of Methanosarcinaceae. Genes horizontally transferred from bacteria to archaea neither tend to be transferred with their regulators, nor in long operons.

CleavPredict: A Platform for Reasoning about Matrix Metalloproteinases Proteolytic Events.

CleavPredict (http://cleavpredict.sanfordburnham.org) is a Web server for substrate cleavage prediction for matrix metalloproteinases (MMPs). It is intended as a computational platform aiding the scientific community in reasoning about proteolytic events. CleavPredict offers in silico prediction of cleavage sites specific for 11 human MMPs. The prediction method employs the MMP specific position weight matrices (PWMs) derived from statistical analysis of high-throughput phage display experimental results. To augment the substrate cleavage prediction process, CleavPredict provides information about the structural features of potential cleavage sites that influence proteolysis. These include: secondary structure, disordered regions, transmembrane domains, and solvent accessibility. The server also provides information about subcellular location, co-localization, and co-expression of proteinase and potential substrates, along with experimentally determined positions of single nucleotide polymorphism (SNP), and posttranslational modification (PTM) sites in substrates. All this information will provide the user with perspectives in reasoning about proteolytic events. CleavPredict is freely accessible, and there is no login required.

Mycobacterial nicotinate mononucleotide adenylyltransferase: structure, mechanism, and implications for drug discovery.

Nicotinate mononucleotide adenylyltransferase NadD is an essential enzyme in the biosynthesis of the NAD cofactor, which has been implicated as a target for developing new antimycobacterial therapies. Here we report the crystal structure of Mycobacterium tuberculosis NadD (MtNadD) at a resolution of 2.4 Å. A remarkable new feature of the MtNadD structure, compared with other members of this enzyme family, is a 310 helix that locks the active site in an over-closed conformation. As a result, MtNadD is rendered inactive as it is topologically incompatible with substrate binding and catalysis. Directed mutagenesis was also used to further dissect the structural elements that contribute to the interactions of the two MtNadD substrates, i.e. ATP and nicotinic acid mononucleotide (NaMN). For inhibitory profiling of partially active mutants and wild type MtNadD, we used a small molecule inhibitor of MtNadD with moderate affinity (Ki ∼ 25 µM) and antimycobacterial activity (MIC80) ∼ 40-80 µM). This analysis revealed interferences with some of the residues in the NaMN binding subsite consistent with the competitive inhibition observed for the NaMN substrate (but not ATP). A detailed steady-state kinetic analysis of MtNadD suggests that ATP must first bind to allow efficient NaMN binding and catalysis. This sequential mechanism is consistent with the requirement of transition to catalytically competent (open) conformation hypothesized from structural modeling. A possible physiological significance of this mechanism is to enable the down-regulation of NAD synthesis under ATP-limiting dormancy conditions. These findings point to a possible new strategy for designing inhibitors that lock the enzyme in the inactive over-closed conformation.

Sequence and conformational specificity in substrate recognition: several human Kunitz protease inhibitor domains are specific substrates of mesotrypsin.

Mesotrypsin is an isoform of trypsin that is uniquely resistant to polypeptide trypsin inhibitors and can cleave some inhibitors rapidly. Previous studies have shown that the amyloid precursor protein Kunitz protease inhibitor domain (APPI) is a specific substrate of mesotrypsin and that stabilization of the APPI cleavage site in a canonical conformation contributes to recognition by mesotrypsin. We hypothesized that other proteins possessing potential cleavage sites stabilized in a similar conformation might also be mesotrypsin substrates. Here we evaluated a series of candidate substrates, including human Kunitz protease inhibitor domains from amyloid precursor-like protein 2 (APLP2), bikunin, hepatocyte growth factor activator inhibitor type 2 (HAI2), tissue factor pathway inhibitor-1 (TFPI1), and tissue factor pathway inhibitor-2 (TFPI2), as well as E-selectin, an unrelated protein possessing a potential cleavage site displaying canonical conformation. We find that Kunitz domains within APLP2, bikunin, and HAI2 are cleaved by mesotrypsin with kinetic profiles of specific substrates. TFPI1 and TFPI2 Kunitz domains are cleaved less efficiently by mesotrypsin, and E-selectin is not cleaved at the anticipated site. Cocrystal structures of mesotrypsin with HAI2 and bikunin Kunitz domains reveal the mode of mesotrypsin interaction with its canonical substrates. Our data suggest that major determinants of mesotrypsin substrate specificity include sequence preferences at the P1 and P'2 positions along with conformational stabilization of the cleavage site in the canonical conformation. Mesotrypsin up-regulation has been implicated previously in cancer progression, and proteolytic clearance of Kunitz protease inhibitors offers potential mechanisms by which mesotrypsin may mediate pathological effects in cancer.

Basis for substrate recognition and distinction by matrix metalloproteinases.

Genomic sequencing and structural genomics produced a vast amount of sequence and structural data, creating an opportunity for structure-function analysis in silico [Radivojac P, et al. (2013) Nat Methods 10(3):221-227]. Unfortunately, only a few large experimental datasets exist to serve as benchmarks for function-related predictions. Furthermore, currently there are no reliable means to predict the extent of functional similarity among proteins. Here, we quantify structure-function relationships among three phylogenetic branches of the matrix metalloproteinase (MMP) family by comparing their cleavage efficiencies toward an extended set of phage peptide substrates that were selected from ∼ 64 million peptide sequences (i.e., a large unbiased representation of substrate space). The observed second-order rate constants [k(obs)] across the substrate space provide a distance measure of functional similarity among the MMPs. These functional distances directly correlate with MMP phylogenetic distance. There is also a remarkable and near-perfect correlation between the MMP substrate preference and sequence identity of 50-57 discontinuous residues surrounding the catalytic groove. We conclude that these residues represent the specificity-determining positions (SDPs) that allowed for the expansion of MMP proteolytic function during evolution. A transmutation of only a few selected SDPs proximal to the bound substrate peptide, and contributing the most to selectivity among the MMPs, is sufficient to enact a global change in the substrate preference of one MMP to that of another, indicating the potential for the rational and focused redesign of cleavage specificity in MMPs.

Characterization of bacterial NMN deamidase as a Ser/Lys hydrolase expands diversity of serine amidohydrolases.

NMN deamidase (PncC) is a bacterial enzyme involved in NAD biosynthesis. We have previously demonstrated that PncC is structurally distinct from other known amidohydrolases. Here, we extended PncC characterization by mutating all potential catalytic residues and assessing their individual roles in catalysis through kinetic analyses. Inspection of these residues' spatial arrangement in the active site, allowed us to conclude that PncC is a serine-amidohydrolase, employing a Ser/Lys dyad for catalysis. Analysis of the PncC structure in complex with a modeled NMN substrate supported our conclusion, and enabled us to propose the catalytic mechanism.